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Abstract — In order to make full use of geographic routing 
techniques developed for sensor networks, nodes must be local- 
ized. However, traditional localization and virtual localization 
techniques are dependent either on expensive and sometimes 
unavailable hardware (e.g. GPS) or on sophisticated localization 
calculus (e.g. triangulation) which are both error-prone and with 
a costly overhead. 

Instead of actually localizing nodes in the physical two- 
dimensional Euclidean space, we use directly the raw distance to 
a set of anchors to produce multi-dimensional coordinates. We 
prove that the image of the physical two-dimensional Euclidean 
space is a two-dimensional surface, and we show that it is possible 
to adapt geographic routing strategies on this surface, simply, 
efficiently and successfully. 

I. Introduction 

If in wired networks each node is equipped with substantial 
computation and storage resources, this is not the case for 
sensor networks which are made of small and cheap devices 
and therefore can not maintain routing tables. Instead of using 
routing tables, local routing techniques have been developed. 
A compelling technique - geographic routing - consists in 
using nodes' coordinates. Many algorithms have been devised, 
such as GPSR |5| and OAFR J6) which use a combination 
of greedy routing and perimeter routing. One can also cite 
GRIC ifTUl a greedy routing algorithm following the sides of an 
obstacle when one is met, and which introduces some inertia in 
the direction followed by the message. Early obstacle detection 
algorithms, that use of a bit of memory at each node, have been 
proposed in J4|, 

Unfortunately, obtaining physical coordinates is problematic 
in itself. The hypothesis of having a GPS for each sensor 
arguably leads to too expensive and heavy devices. This 
assumption may be weakened by equipping some nodes called 
anchors with the GPS or with other localization hardware. 
Approximate coordinates are then computed for all nodes 
of the network in a localization phase. In Q, three such 
algorithms are compared, namely Ad-hoc positioning, Robust 
positioning, and N-hop multilateration. One can also cite 
the algorithm At-Dist [13], which is a distributed algorithm 
estimating the position of each node together with an estimate 
of its accuracy. Some authors improved these results by using 
angle measurements [1]. All of these localization techniques 
invariably require a flooding from the anchors and many com- 
putations at each node. They are therefore energy consuming, 
error-prone and compute only approximate coordinates. 

Another approach followed in Q, (8), ifTTl . Ifl2l consists 
in computing virtual coordinates and has the advantage of not 
needing anchors. In IfTTl . the authors use a mobile unit to 
assist in measuring the distance between nodes in order to 
improve accuracy. The algorithm proposed in 12) first chooses 



three nodes that will behave as anchors and from which virtual 
coordinates will then be determined. If these techniques do 
not need any external hardware, they also suffer from their 
inaccuracy or high energy consumption in a preprocessing 
phase. 

In this paper, we discard any preprocessing technique and 
propose to directly use raw distance information. We study 
routing algorithms using directly the distance to the anchors as 
coordinates, as first proposed in (3), without computing from 
them 2-dimensional coordinates. In section [II] we precisely 
describe how the idea is implemented, and in section [III] we 
analyze how a message sent towards a destination performs in 
the new coordinate system. 

II. Implementation 

Current localization methods rely on raw information com- 
puted externally from normal sensing nodes (exact location 
of some anchors), and on raw information computed locally 
in normal sensing nodes (distance to anchors, angle mea- 
surements). In this paper, we do use the information about 
the distance to some anchors, but we completely discard 
any physical information that the anchors might have. This 
gives much more flexibility in the way sensor networks are 
deployed: anchors might be some external entities, as a plane 
or a robot; anchors might be specialized nodes whose only 
purpose is to emit a strong signal, or they might be randomly 
chosen sensors which advertise their distance to the other 
nodes. 

We build a multi-dimensional coordinate system using di- 
rectly the raw information, i.e. the distance to the anchors. 
Given a node at location X, we define the multi-dimensional 
coordinates f(X) of this node as its distance to the anchors 
at location Ai, A 2 , . . . A n : 



(d{X,A x )\ 
d{X, A 2 ) 

\d(X,A n )J 



We call this function the anchor coordinates function, and 
we call these multi-dimensional coordinates the anchor coor- 
dinates. Whereas any distance function, such as hop count, 
may be used [3|, in section III we pay a special attention to 



the properties of / when d is the Euclidean distance. 

In the next subsection we discuss the computation costs that 
are specific to using multi-dimensional coordinates. We then 
go into the details of greedy routing implementation, and into 
the details of rotating multi-dimensional vectors. 



A. Computation Cost 

While saving on initialization overhead, multi-dimensional 
routing causes some additional computation costs when send- 
ing messages compared to traditional two-dimensional routing. 
Here is a break-down of various vector operations: 
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Note that additions and multiplications typically use 1 CPU 
cycle, whereas the expensive operations (square root extrac- 
tion, inversion) stay the same in multi-dimensional routing 
as in traditional two-dimensional routing. We also point out 
that theses computation costs are not communication costs and 
are lower in terms of energy consumption by some order of 
magnitude. 

B. Greedy Routing 

Greedy routing is the most basic geographic routing algo- 
rithm. It consists in following the direction to the destination. 
This basic strategy is widely used as a default mode in most 
geographic routing protocols. When a node at location X 
which wants to send a message towards a final destination 
at location D, three implementations of greedy routing are 
routinely used: 

1) (canonical) for each neighbor location X', compute the 

distance d(X' , D) and send the message to the neighbor 

— - — y — - — y 

which is closest to D. Alternatively, compute X D-X D 
instead of d(X',D). 

2) for each ne ighb or location X', compute the scalar prod- 
uct XX' ■ X~t> and select the neighbor with the best 
result. 

3) for e ach n eighbor location X', compute the scalar prod- 



uct 

\\xx'\\ 
result. 



xf> and select the neighbor with the best 



These three implementations are valid for any number of 
coordinates. 

C. Rotation 

When greedy strategies fail, a number of two-dimensional 
routing algorithms fall back on more sophisticated routing 
modes that use rotations or angle computations Q, |10|. 
When using two dimensions, a rotation is typically defined by 
rot a : (x, y) — > (x cos a + y sin a, y cos a — x sin a). We can't 
define such a rotation in n dimensions (n > 3). However, if we 
assume that our sensors were on a two-dimensional physical 
plane in the first place, then they are distributed over a two- 
dimensional surface in the multi-dimensional space (more on 
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1) compute an orthonormal basis 
plane in f(X) (see section [5T1. 

2) express vectors U as x u i +_^« j + e u by computing 
x u = it ■ i and y u = u ■ j . We assume that it is 
close to the tangent plane in f(X), which means that 
we ignore in fact e^. 

Rotations are then normally carried out on the tangent plane. 
The sensitive part is to compute ( i , j ) and to make sure that 
the orientation of the surface is preserved when routing the 
message (taking the surface upside-down has the undesirable 
effect of negating angles). Given a node at location X, a 
destination at D, and a basis (i id, joid) inherited from a 
previous node, we do the following: 

1) choose two neighbors at position X\ and X2 

• either arbitraril y (low qu ality, inexpensive) 

I XX ^ XX ^ I 

• or such that - — is minimal (i.e. choose 

\jxxl\\\\xxl\\ 



2) compute i 

3) compute it 



XX\ and XX2 as orthogonal as possible) 
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4) compute ~3 = -\^r\ 

5) compute a = (i ■ ( oW ) (it ■ j old ) 
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6) if a > then set j = it, else set j = —it 
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Note that many algorithms using angles use normalized 
vectors. Therefore,_mqst of the normalization cost when com- 
puting the basis ( i , j ) is not an additional cost of multi- 
dimensional routing. 

D. Experimentation 

The routing protocols GRIC OH and ROAM gl were 
experimented in J3|, using 6 anchors randomly distributed in 
the network. Messages were delivered with 100% success, 
and the average length of the paths followed by messages 
using anchor coordinates was the same as when using physical 
Euclidean coordinates. 

III. Algebraic Analysis 

In the plane with Euclidean distance, any node has a 
pair of physical coordinates X = {x,y). We denote by 
Ai = (xi,yi) the physical coordinates of the i th anchor. The 
anchor coordinates function is a function from R 2 — > K n 
defined by 



/ ^(x - Xl ) 2 + (y - Vl ) 2 \ 
y/(x - x 2 ) 2 + {y - y 2 ) 2 

V^/ (x - x n ) 2 + (y- y n ) 2 ) 



this in section IIIi. We do the following: 



Since the functions fi : (x,y) — > yj(x — Xi) 2 + (y — y^) 2 are 
continuous and C°° except in (xi, yt), we show that with three 
or more anchors that are not on the same line, the image /(M 2 ) 
in E" is a continuous surface (claim [TJ. Figure [TJ represents 
the image of /, when there are three anchors at location (0,0), 
(0,1) and (1,0). 



Surf ace uith 3 anchors 
sqrt(u**2+<v-l>**2), sqrt(u**2+v**2>, -sqrt((u-l>**2+y**2> 
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Fig. 1. Representation of the distance to three anchors. 



First, we describe in subsection III-A the vector spaces 
that are tangent to /(M 2 ). Next, we express in subsection 



III-B what is the physical direction of messages that use the 
greedy strategy with virtual coordinates. This physical direc- 
tion produces a curve that approximates the paths followed 



by messages. We then tell in subsection III-C what are the 
convergence conditions on /(M 2 ) under which the curve ends 
at the destination, and prove a bound on the length of this 



curve. Finally, we study in subsection III-D how the placement 
of anchors affect the convergence conditions and how we can 
guarantee that they are met. 

A. Tangent Space 

At any point f(x, y), the surface /(M 2 ) has a tangent vector 



space spanned by the two vectors ^(x,y) and ?f-(x, y). We 
have 
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Claim 1 The vector space that is tangent to the surface /(R 2 ) 
in f(X) is two-dimensional if and only if the node X and the 
anchors A\, A2, . . . ,A n are not situated on a single line in 
the physical space. 



Proof: The tangent vector space is two-dimensional if and 
only if fa(x,y) and ^(x,y) are not collinear. Conversely 
^(x,y) and g(x, y) are collinear if and only if there is 
a € [0, 2ir[ such that ^(x,y) cos a + y) sin a = 0. By 

changing the physical coordinates into u = x cos a + y sin a 
and v — y cos a — x sin a (we also set Ui = Xi cos a + yi sin a 
and Vi — yi cos a — Xi sin a), we express the tangent vector 
space with the two vectors 
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V — v n 
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We have |^ (X) = if and only if for all i G {1 . . . n}, u 



When fa:(x,y) and ^-(x, y) are not collinear, then the 
Jacobian matrix 
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\ \J (x~x n ) 2 + (y-y n ) 2 yj (x-x n ) 2 + (y-y n ) 2 

defines a morphism of the physical plane into the vector space 
tangent to f(M. 2 ) at f(x, y). Given a node at position X in the 
physical space and its neighbors at position X±,X2, ■ ■ ■ ,Xg, 
it is not unreasonable to assume that for all i, f(Xi) is close 
to the Taylor expansion f(X) + J j{X){XX > i ) in the affine 
space tangent to /(M 2 ) in f(X). 

B. Directional Vector 

In a greedy routing strategy using virtual coordinates, the 
neighbor X' of choice will be a maximum for some scalar 



product f(X)f(X') ■ f(X)f(D). 

Claim 2 Given two physical positions X, D € M 2 , the 
function sx : M 2 — > M such that for any vecto r XX ' E 
R 2 , sx(X X') is the scalar product of J f (X)(XX') by 
f(X) f(D) is a linear form that can be expressed as 



XX' -> XX' ■ a iXAi 



where a,, 



,Ai)-d(D,Ai) 
d(X.Ai) 



Proof: The transformation XX' -)• J f (X)(XX') is a 
linear function. Since the scalar product by f(X)f(D) is a 
linear form, sx is also a linear form. We may decompose the 



vector f(X)f(D) into £\ (d{D, AA - d(X, A,))!, where 1, 
is the multi-dimensional vector with 1 as its i th coordinate and 
zeroes everywhere else. In this manner, sx — Yli s x,i where 



8 X ,i{XX') = {d{D, A,) - d(X, A^JfiX^XX 1 ) ■ U 



J f (X)(XX') ■ u = 



(x - Xi)(x' -x) + (y- yi)(y' - y) 



^J{x- xA 2 + (y - j/j) 2 
Thus sx.i can be expressed as 
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Fig. 2. Contribution of anchors in Pi, P2, P3, P4. 
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Given a node at physical location X and a destination D G 
I 2 , we call apparent destination related to D in X the location 



D' 



X + Y^ Q!iX% = X + Y 



d(X,A i )-d(D,A i )—> 



d{X, Ai) 



XAi 



C. Virtual consistency 

We say that the anchor coordinate system is virtually 
consistent at distance r for a physical destination D G M 2 , 
if at every point I/fl such that f{X) is in a closed metric 
ball of center f(D) and radius r, then sx 7^ 0. Note that 
s x = if and only if the multi-dimensional vector f(X)f{D l ) 
is orthogonal to the vector space tangent to f(M. 2 ) in f(X). 
It is also equivalent to state that the anchor coordinate system 
is virtually consistent at distance r for a physical destination 
D G W 2 , if no closed metric ball centered on f(D) and of 
radius < r' < r is tangent to /(M 2 ). 

Claim 3 If the anchor coordinate system is virtually consis- 
tent at distance r for a physical destination D € IR 2 , then 
there is X G R + such that for any point Xq with f(Xo) in 
a closed metric ball of center f(D) and radius r we have a 
curve c[0, 1] G R 2 that verifies: 

• c : [0, 1] — > M. 2 is a derivable function, 
. c(0) = X and c(l) = D, 

• At any point t G [0, 1[, the vector ||| (t) is collinear with 

the vector c(t)D[ where D' t is the apparent destination 
related to D in c(t). 

• fS\\^\\dt<Xd(X ,D). 

Proof Let k be the largest positive number such that 
for any point X = (x,y) with f(X) in a closed metric ball 
of center f(D) and radius r, the orthogonal projection of 



f(X)f(D) on the vector space defined by the two vectors 
gjOc, y) and ^(x,y) has a norm greater than or equal to 
kd(X, D). Since / is a continuous function, the set of physical 
positions X such that d(f(X),f(D)) < r is compact subset 
of M 2 . Therefore, if k was equal to z ero^ then there would be 
a point X ^ D in the ball such that XL) is orthogonal to the 
surface f(M. 2 ), which we excluded in our assumptions. 

Let c : [0, 1] -> M 2 be the function defined by c(0) = 
Xq and such that g° c ^ (t ) is the orthogonal projection of 



k d((fo f cm,/( i D)) ) (/ c )( t )/( I) j on the vectOT s P ace defined 
by the two vectors §£(c(i)) and |^(c(t)). Since 
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we can see that 
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(t)>d(f(X ),f(D)) 



which implies that c(l) = D. The norm of ff (i) is smaller 
than or equal to ||(J / (c(<)))- 1 ||d(/(X ), /(£>)), which means 
that 



|J^|M*< max MJMtyr^W&oJiD))) 



I — ||dt < max IKJ^cW))- 1 !!^,^). 
<9(*) «e[o,i] 



D. Physical consistency 

We say that the anchor coordinate system is phy si cally con- 
sistent at position X for the destination D if XD' ■ XD > 0, 
where D' is the apparent destination related to D' in X. 
Observe that if the anchor coordinate system is physically 
consistent for the destination I? in a ball B around D, then it 
is virtually consistent at distance r for the physical destination 
D, where r is the radius of the biggest multi-dimensional ball 
n such that Un/(R 2 )c f{B). 

To study the physical consistency of the system at position 
X for the destination D, we split the phys ical plane in 
four parts Pi,P9,P 3 ,F4 with Pi = {X'\XX' -XD<0}, 
P 2 = {X'\XX' ■ X~f) > and d(X,X') < d(X,D)}, 
P 3 = {X'\DX' ■ Dlt > and d(X,X') > d(X,D)}, 
< 0}. Since the apparent destination 

D' is defined by 



D' = X 



d(X,A i )-d(D,A i 
d{X,AA 



XAi 



we see as illustrated in Figure [2| that only the anchors in P2 
give a negative contribution to XD' ■ XD. 



If anchors are randomly distributed in the network, the neg- 
ative contribution will most probably be small enough for the 
system to be consistent, unless Pi and P4 are almost void of 
nodes, which happens when X and D are located on opposite 
borders of the network (so that all the anchors are between 
them). This situation never occurred in the experiments carried 
out in 1 3 1 (but the end nodes were not chosen on the border 
of the network). Nevertheless, physical inconsistency may be 
avoided by selecting anchors when the destination D of a 
message originating from Xq is far away: 

1) by default, use all the anchors. 

2) compute^ = maxj 6 { lj max(d(D, Ai), d(Xo, Ai)) = 
max(||/(L>)|| 00 , ||/(Xo)||oo)- U gives an idea of the 
diameter of the network. 

3) for each node X along the path of the message 

a) compute l x = maxj e /i n \\d(X,D) — 
d(X,Ai)\ = \\f(D) - f(X)\Q l x is smaller 
than d(X, D). 

b) if using all the anchors and if lx > ^ then use 
only the anchors Ai such that d(D,Ai) < 

c) if using a subset of anchors and if lx < then 
use all the anchors. 

In this way, physical inconsistency can be completely avoided 
in the network, at the cost of using a different coordinate 
system when d(X, D) is comparable to the diameter of the 
network. 
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IV. Conclusion 



Geographic routing is an essential component in connecting 
sensor networks. Foregoing the previously necessary localiza- 
tion phase where physical Cartesian coordinates are produced 
is an important step into making networks more robust and to- 
tally independent from external hardware. Sensor network ap- 
plications that use localization information exclusively inside 
the network may transparently use virtual coordinates, whereas 
sophisticated physical localization may still be performed 
at some external base station from the virtual coordinates 
whenever localization must be used externally. In this way, 
directly using raw distance information without any costly 
or sophisticated localization calculus is a simple, viable, and 
efficient way to perform geographic routing. 
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